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Galaxies congregate in clusters and along filaments, and are miss- 
ing from large regions referred to as voids. These structures are seen 
in maps derived from spectroscopic surveys’” that reveal networks 
of structure that are interconnected with no clear boundaries. Ex- 
tended regions with a high concentration of galaxies are called ‘super- 
clusters’, although this term is not precise. There is, however, another 
way to analyse the structure. If the distance to each galaxy from Earth 
is directly measured, then the peculiar velocity can be derived from 
the subtraction of the mean cosmic expansion, the product of distance 
times the Hubble constant, from observed velocity. The peculiar ve- 
locity is the line-of-sight departure from the cosmic expansion and 
arises from gravitational perturbations; a map of peculiar velocities 
can be translated into a map of the distribution of matter’. Here we 
report a map of structure made using a catalogue of peculiar velo- 
cities. We find locations where peculiar velocity flows diverge, as 
water does at watershed divides, and we trace the surface of diver- 
gent points that surrounds us. Within the volume enclosed by this 
surface, the motions of galaxies are inward after removal of the mean 
cosmic expansion and long range flows. We define a supercluster to 
be the volume within such a surface, and so we are defining the extent 
of our home supercluster, which we call Laniakea. 

The distribution of matter can be determined by two independent 
methods: either based on surveys of the distribution of galaxies in pro- 
jection and redshift, or from the motions of galaxies. With the former, 
using galaxy redshift surveys, the assumption is required that the galaxy 
‘lighthouses’ and mass distribution are strongly correlated, a condition 
that requires confirmation if, as is suspected, only a minor fraction of 
matter is baryonic. Moreover, with the former there is a stringent demand 
that the survey be complete, or at least that its incompleteness be well 
understood. With the latter, studies of galaxy motions, sparse sampling 
is acceptable (indeed inevitable) but dealing with errors is a challenge. 
Except for the very closest galaxies, uncertainties in distance measure- 
ments translate into uncertainties in the peculiar velocities of galaxies 
that are larger in amplitude than the actual peculiar velocities. Many mea- 
surements are required for suitable averaging and care must be taken to 
avoid systematic errors. Overall, the two paths to determining the dis- 
tribution of matter are in good agreement, a consequence that represents 
a considerable success for the standard model of structure formation 
via gravitational instability“. 

The path from velocities to mass distributions benefits from the coher- 
ence in velocities on large scales. Multipole components in the velocity 
field can point to tidal influences beyond the survey region. The current 
all-sky redshift surveys and distance measurement surveys reach sim- 
ilar depths, but the latter probe structure to greater distances because of 
sensitivity to uncharted attractors and repellers. Coherence in motions 
on large scales means that signals can be measured by averaging over 
data in circumstances where individual contributions are very noisy. 

Details about the actual measurement of galaxy distances and the der- 
ivation of peculiar velocities are given in Methods. These two parameters 
are available for more than 8,000 galaxies, affording extremely detailed 
information locally, degrading outward to increasingly coarse coverage. 
We use this material to reconstruct the large-scale structure of the near- 
by Universe’. 


We obtain the underlying three-dimensional velocity and density fields 
by the Wiener filter algorithm*’, assuming the standard model of cos- 
mology as a Bayesian prior. Large-scale structure is assumed to develop 
from gravitational instabilities out of primordial random Gaussian fluc- 
tuations. The developing density and velocity fields retain their Gauss- 
ian properties as long as the growth is in the linear regime. It has been 
shown that with a random Gaussian field, the optimal Bayesian estima- 
tor of the field given the data is the Wiener filter minimal variance esti- 
mator. At the present epoch, large-scale structure has become nonlinear 
on small scales. However, it is an attractive feature of the velocity field 
that the break from linearity is only on scales of a few megaparsecs, an 
order of magnitude smaller in scale than the deviations from linearity 
for the density field. In any event, the present discussion concerns struc- 
ture on scales of tens to hundreds of megaparsecs, comfortably in the 
linear regime. 

The Wiener filter result is determined by the ratio of power to power 
+noise. Hence, the large-scale structure is strongly constrained nearby, 
where uncertainties are small and the coverage is extensive. At large 
distances, where the data become more sparse and noisy, the Wiener 
filter attenuates the recovered density and velocity fields to the null field 
that is expected in the absence of data. However in the linear regime 
there is coherence in galaxy flows on much larger scales than seen in 
density fluctuations. Tidal influences from beyond the surveyed regions 
can be manifested in cosmic flows on scales that exceed the coverage in 
measured distances by a factor of two (ref. 10). 

The ultimate goal is to map the velocity field to a radius that com- 
pletely encompasses the sources of the motion of the Local Group (of 
631 km s`’) reflected in the cosmic microwave background dipole". How- 
ever, our knowledge of flows on large scales remains inevitably modulated 
by the extent of the data. Our analysis of the data in the Cosmicflows-2 
galaxy catalogue (see Methods for details) tells us that a coherent flow 
extends across the full extent of the region that we can map, reaching 
the Shapley concentration”. It is clear that we do not yet have a suffi- 
ciently extensive compendium of distances to bound the full source of 
our deviant motion from the cosmic expansion. 

For the present discussion, we focus on intermediate-scale flow pat- 
terns. The standard model of cosmology predicts that on the scale we 
are considering the flow is irrotational: namely, the velocity field v is 
the gradient of a potential ¢, v = — V@. The local minima and maxima 
of the potential (attractors and repellers respectively) are the drivers of 
the large-scale flow. We can define a ‘basin of attraction’ as the volume 
containing all points whose flow lines converge at a given attractor. The 
large-scale structure can be characterized on scales of a few megapar- 
secs and above by attractors and their basins of attraction. 

The Wiener filter provides a straightforward way of decomposing the 
velocity field into a local component that is induced by the distribution 
within a zone and a tidal residual". In the linear regime, the velocity 
and density fields, vand 6, are directly related: V- v= — Hof(Qm; 24) 
where f depend on the cosmological matter and vacuum energy densi- 
ties characterized by Qm and Q 4. We specify a centre and radius, and 
the density field within the defined volume is assumed to be a cut-out 
from the full Wiener filter density map. Then the Poisson-like equa- 
tion between velocity and density can be solved to derive the velocity 
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field responding to just the matter within the prescribed volume. The 
vector subtraction of the local velocity component from the full flow 
gives the external component of the velocity field. The residual com- 
ponent is responsible for the bulk motion of the zone under considera- 
tion and for a quadrupole component within the zone. The decomposition 
allows us to probe the local velocity field, with the tidal field induced 
by distant structures filtered out. Relative attractors and their basins of 
attraction are defined with respect to that local field. 

One more useful tool is now mentioned before we turn to results. At 
each position in space, the three eigenvalues of the velocity shear tensor 
can be calculated. If these eigenvalues are ordered from most positive 
to most negative, then a threshold can be set that captures four possi- 
bilities. Flows can be inward on all three axes, the condition ofa cluster, 
inward on two axes and outward on the third, the condition of a fila- 
ment, inward on one axis and outward on two, hence a sheet, or out- 
ward on all three axes, hence a void. Boundaries can be created around 
contiguous regions with the same shear properties, and the contours out- 
line the cosmic web as reconstructed by the V-web algorithm‘. 

We now consider the full presentation of our results given in the Sup- 
plementary Video and in Extended Data Figs 1-5. We pay attention in 
turn to smaller (but still in the linear regime) scales to examine the sep- 
aration of local and tidal flows and to isolate local basins of attraction. 
Particular attention has been given to locations where there are local di- 
vergences. In co-moving coordinates and with the removal of long-range 
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Figure 1 | Two views of the 
Laniakea supercluster. The outer 
surface (blue) demarcates the limits 
of local velocity flows. The plot is in 
supergalactic coordinates with axes 
SGX and SGY shown in the plane at 
SGZ = 0. Our Milky Way galaxy is at 
the origin. Units of 1,000 kms™' in 
velocity correspond to roughly 

13.3 Mpc. Velocity streamlines are 
shown in black and terminate in the 
vicinity of the Norma cluster. 
Individual galaxies from a redshift 
catalogue are given colours to 
distinguish major components 
within the Laniakea supercluster: the 
historical Local supercluster in green, 
the Great Attractor region in orange, 
the Pavo-Indus filament in purple, 
and structures including the Antlia 
wall and Fornax-Eridanus cloud in 
magenta. Several major entities are 
named. Norma, Hydra, Centaurus, 
Virgo, Ophiuchus, A2870, A3581 
and A3656 are individual clusters of 
galaxies embedded within the 
Laniakea supercluster. Shapley, 
Hercules, Coma and Perseus—Pisces 
are complexes of galaxies outside 
Laniakea. The outer black circle 
defines the domain used to separate 
between local and tidal flows. The 
panels provide two perspectives of 
the same scene. 


l 
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flows, there are places where relatively neighbouring galaxies can be 
found to be moving in opposing directions towards separate local basins 
of attraction. Voids are usually the demarcations between attraction 
basins, but divergences can occur along filaments and sheets. We note 
the very nearby case between our home basin of attraction and the 
Perseus—Pisces complex’>"®, where we find a particularly impressive ex- 
ample of velocity flows in an apparent bridge between attractors that 
are diverging. This divergence occurs in the feature called the Arch in 
Extended Data Fig 2. Similar structures abound on close inspection. 
Velocity information reveals the locations of divergence along filaments 
between high density regions. This dissipation of the cosmic web is ex- 
pected with the accelerated expansion of the Universe. We emphasize 
that peculiar velocity information can reveal details that are otherwise 
hard to discern. 

The particular interest in the present discussion is with the largest 
structure that can be circumscribed within the currently available dis- 
tance and peculiar-velocity data: the structure schematically illustrated 
in Figs 1 and 2. The region includes 13 Abell clusters (with the Virgo 
cluster). Local flows within the region converge towards the Norma and 
Centaurus clusters, in good approximation to the location of what has 
been called the ‘Great Attractor’’”. This volume includes the historical 
Local and Southern superclusters’’, the important Pavo-Indus filament, 
an extension to the Ophiuchus cluster, the Local Void, and the Sculptor 
and other bounding voids. This region of inflow towards a local basin 
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of attraction can be reasonably called a supercluster. The region, if 
approximated as round, has a diameter of 12,000 km s~ ' in units of the 
cosmic expansion or 160 megaparsecs, and encompasses ~10'” solar 
masses. We propose that this region be named the Laniakea super- 
cluster of galaxies (from the Hawaiian; lani, heaven, and akea, spacious, 
immeasurable). 


Online Content Methods, along with any additional Extended Data display items 
and Source Data, are available in the online version of the paper; references unique 
to these sections appear only in the online paper. 
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Figure 2 | A slice of the Laniakea 
supercluster in the supergalactic 
equatorial plane. Shaded contours 
represent density values within the 
equatorial slice, with red at high 
densities, green at intermediate 
densities and blue in voids. Our 
Milky Way galaxy is located at the 
black dot at the origin of the 
supergalactic coordinates system: a 
red arrow points right from the 
black dot toward increasing SGX 
and a green arrow points up toward 
increasing SGY. Individual galaxies 
from a redshift catalogue are given 
as white dots. Velocity flow 
streams within the Laniakea basin 
of attraction are shown in white. 
The orange contour encloses the 
outer limits of these streams. 

This domain has a extent of 
~12,000 km s7 }(~160 Mpc 
diameter) and encloses ~10"7 
solar masses, Mo. 
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METHODS 


The present discussion draws on a new catalogue of galaxy distances and peculiar 
velocities’, one that extends to recession velocities of 30,000 kms’! (redshift z = 
0.1) and with 8,161 entries provides a density of coverage previously unknown. The 
new catalogue is called Cosmicflows-2 and the six main methodologies for the dis- 
tance estimates rely on the characteristics of Cepheid star pulsations, the luminosity 
terminus of stars at the tip of the red giant branch, surface brightness fluctuations 
of the ensemble of stars in elliptical galaxies, the standard candle nature of super- 
novae of type Ia, the adherence by elliptical galaxies to a fundamental plane in lumi- 
nosity, radius, and velocity dispersion, and the correlation between the luminosities 
of spirals and their rates of rotation. Each of the methodologies has strengths and 
weaknesses. The Cepheid and tip of the red giant branch techniques provide high 
precision distances but only to very nearby galaxies. The elliptical fundamental 
plane and spiral luminosity—rotation methods provide individually less accurate 
distances but can be used to acquire samples of thousands of galaxies in a large 
volume. Type Ia supernovae are excellent distance indicators and can be seen far 
away but they arise serendipitously and current samples are small. Jointly the sky 
coverage is now substantial within ~100 Mpc, with only spotty coverage out to 
~400 Mpc. 

When the Cosmicflows-2 compendium was being compiled, considerable effort 
went into ensuring that the six independent methodologies were on acommon scale. 
This scale determines the value of the Hubble constant. In Cosmicflows-2, this pa- 
rameter is determined from the velocities and distances of type Ia supernovae at 
redshifts 0.03 < z <0.5, ina domain where peculiar velocities should be a negligible 
fraction of observed velocities. We found”? Hy = 75.2 + 3.0kms ‘Mpc 1. 

The value of Ho remains somewhat contentious, especially with the value of 67.3 
+ 1.2kms~ | Mpc’ claimed by the Planck collaboration”!. Other determinations 
that actually measure distances and velocities are more compatible with our value””’. 
Variations of the assumed Ho within the Cosmicflows-2 analysis introduce a mono- 
pole term: infall if Ho is increased and outflow if Ho is decreased. Our value mini- 
mizes the monopole term. Tests with the current distances demonstrate unphysical 
monopole terms if Ho is varied by more than +1.5kms ' Mpc’ from the fiducial 
value. The Planck determination at face value introduces a massive outflow. It is 
to be emphasized that, for a study of galaxy flows such as ours, the value assumed 
for Ho needs to be internally consistent with the scaling of distance measures. A 
rescaling, say by a re-evaluation of the distances to the nearest galaxies, leaves pe- 
culiar velocities unchanged. The extremely low Ho value found by the Planck col- 
laboration is only plausible in the face of the monopole implication if distances to 
nearby calibrator galaxies are significantly greater than currently measured. 

The Wiener filter has been used before to reconstruct the underlying flow field 
from sparse and noisy velocity surveys™*”’. Here it is used to overcome the Malmquist 
bias introduced by errors in distance that scatter galaxies out of regions containing 
more galaxies towards emptier regions”. This correction is done in two steps. In 
the first, the Wiener filter is applied to the raw data. New estimated distances are 
defined by dwr = (Vobs — Vwe)/Ho where Vwris the radial component of the Wiener 
filter velocity field at the measured position. This adjustment brings the estimated 
distances to be much closer to their position in redshift space. The Wiener filter is 
then applied to the revised catalogue to yield the reconstructed flow that is presented 
here. The Wiener filter correction of the Malmquist bias has been tested against 
mock catalogues drawn from N-body simulations constrained by the Cosmicflows- 
2 data to reproduce the actual Universe. The Wiener filter strongly suppresses the 
spurious (negative) monopole term introduced by the Malmquist bias. 

The Wiener filter provides the Bayesian mean field, given the data, its uncer- 
tainties and an assumed prior model. The robustness of the Wiener filter recon- 
struction is gauged by sampling the distribution of the residual from the mean field, 
assuming that the underlying field is Gaussian”. To meet this requirement, an en- 
semble of 20 constrained realizations of the underlying flow has been constructed”. 
A local/tidal decomposition has been applied to a sphere of radius 6,000 km s~! 
= 80 Mpc centred on the putative centroid of the attractor region ([—4,700, 1,300, 
—500] kms). Each one of these realizations has been found to harbour a distinct 
monolithic over-dense structure that dominates that sphere. Extending the radius 
of the sphere to a radius of 7,000 km s | and above we find the breakdown of the 
monolithic structure and the emergence of more than one attractor within the 
sphere. The ensemble of constrained realizations can be used to assess the stat- 
istical significance of our claim for the existence of a basin of attraction on a scale 
of 6,000 kms! = 80 Mpc. All of the 20 realizations, each consistent with the data 
and the assumed prior ACDM model, reproduce such an attractor. 

A superiority of the recovery of structure from peculiar velocities with the Wiener 
filter over reconstructions from redshift surveys is the sensitivity to structures in 
regions that cannot be observed, in particular in the heavily obscured equatorial 
band of the plane of our Galaxy. Indeed a large part of our watershed is hidden by 
obscuration. The connections between the important Pavo-Indus and Centaurus 


components of Laniakea, the extension up to Ophiuchus cluster, the filament from 
the Antlia Wall that ultimately connects to the Perseus—Pisces complex, and the 
continuation of the Perseus—Pisces filament beyond Perseus are all examples of hid- 
den structures manifested in flow patterns. The two dominant attractors in Shapley 
and Lepus are both at low galactic latitudes. A further nice example is the identifi- 
cation of the previously unknown Arch. 

The Supplementary Video and the figures presented in this paper were created 
using SDvision”, an interactive visualization software developed within the frame- 
work of IDL Object Graphics. In Supplementary Information we give the transcript 
of the dialogue in versions that include a sound track. The transcript is also pro- 
vided in the closed caption version (transcript of the dialogue embedded in the video). 

The Supplementary Video begins with a half-rotation display of galaxies from 
the 2MASS Extended Source Catalog” with distances given by measured velocities 
assuming the Hubble law. This display gives a general idea of the structure in the 
nearby Universe although there is increasing incompleteness with distance from our 
position at the centre and there is some inhomogeneity in coverage with direction. 

During a second half rotation of the scene, there is a display of the peculiar ve- 
locity vectors that are fundamental to the current analysis. Peculiar motions towards 
us are in blue and those away are in red. Measurements are dense within 7,500 km s j 
= 100 Mpc and fall off quickly beyond that distance. 

The rotations continue, but now with the full velocity flows derived from the 
Wiener filter analysis. The immediate product of the filtering is a most probable 
three-dimensional peculiar velocity at each position within the volume of the study. 
An imaginary vehicle that starts from any seed location is passed to an adjacent po- 
sition by the locally inferred vector of motion, then again passed onward, creating 
flows that are captured by the visualized stream lines. The dominant flow is towards 
the Shapley concentration of galaxies. A second feature involves an accumulation of 
flows in the Perseus-Pisces region that then streams down to the Lepus region. After 
a rotation, the implied large-scale potential field is superimposed. If the viewing 
box is extended, it is reasonably convincing that the flow to Lepus continues on to 
connect with Shapley. In other words, essentially the entire volume being displayed 
is involved in a flow towards Shapley. This proposition will be explored elsewhere. 

The emphasis of the present discussion is on intermediate scales that are within 
the domain of most of the distance measurements. The contours that become su- 
perimposed on the flow lines represent the density field derived from the Wiener 
filter velocity field. The amplitudes of the density reconstructions are high near the 
centre of the cube where data constraints dominate errors, and taper to low values 
at the edges where uncertainties are large and the derived density field approaches 
the mean field. Extended Data Fig. 1 is extracted from this sequence. 

After some more rotations, the flow lines are seen to break up into disjoint pieces. 
These disconnections occur because we have begun to use the trick of separating 
into local and tidal components. Local velocity flows are due to the mass distri- 
bution in a restricted region and tidal flows are due to mass outside the restricted 
region. 

Next there is a transition to the V-web representation (see main text). The eigen- 
values of the velocity shear can be evaluated at each local position. Surfaces are shown 
that enclose clusters, filaments, and sheets. Now in the video there is a zoom inward 
and flow lines are shown for local motions within a zone extending to 6,000 km s~ i 
= 80 Mpc from our central position. We are sampling two separate basins of attrac- 
tion. The red flow lines are directed toward the Perseus-Pisces complex while the 
black flow lines are directed toward the Centaurus-Norma-Pavo-Indus structures. 
Our Earth lies within the domain of the black flow lines, although we are near the 
boundary where local flow directions flip. It is seen that the flows converge to follow 
along filaments. There are interesting places along the V-web filaments where local 
peculiar velocities flip in direction. After some time in this video sequence there is 
the superposition of individual galaxies from the redshift catalogue. Colours are 
given to members of different structures. The positions of the individual galaxies 
are not corrected for redshift distortion but reasonable agreement is seen between 
the locations of galaxies, the flow lines, and the inferred V-web structures. Extended 
Data Fig. 2 is extracted from the video as an illustration of this point. 

In the concluding section of the Supplementary Video, the local/tidal decom- 
position is shifted to a centre near to the Norma cluster in order to isolate our basin 
of attraction. The flow lines now only include those that are moving inward. A sur- 
face is created at the boundary of this region. Figure 1 shows two projections from 
Supplementary Video frames and Fig. 2, the concluding frame of the Supplemen- 
tary Video, is a projection onto the equatorial plane in supergalactic coordinates. 
This latter figure is augmented in Extended Data Fig. 3 to include flow lines away 
from, as well as towards, the local region of attraction. In addition there is an indi- 
cation in Extended Data Fig 3 of the swath of the obscuration in the plane of the 
Milky Way. Extended Data Figs 4 and 5 give two orthogonal views of the inward 
and outward local flows. 

The orange outline in Fig. 2 and in Extended Data Figs 3-5 is the projected bound- 
ary of our local basin of attraction that we call the Laniakea supercluster. This domain 
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of diameter 160 Mpc, given the mean density of the Universe, encloses 10” Mo. 
We stress that in the fullness of distance measurements on a much larger scale it 
will almost certainly be demonstrated that Laniakea is not at rest with respect to 
the cosmic expansion and is only a part of something very large that is at rest in the 
cosmic reference frame. We remind those unfamiliar with the field that cosmic ex- 
pansion velocities are removed in deriving peculiar motions. Infalling motions on 
large scales are only perturbations. All galaxies except those in the immediate vicinity 
of clusters and groups are flying apart. 
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Extended Data Figure 1 | Structure within a cube extending 16,000 kms" * 
(~200 Mpc) on the cardinal axes from our position at the origin. Densities 
ona grid within the volume are determined from a Wiener filter reconstruction 
based on the observed velocity field. Three isodensity contours are shown. 
The density map is detailed near the centre of the box where observational 
constraints are dense and accurate, but tapers to the mean density as constraints 
weaken. Nevertheless, velocity flows illustrated by the black threads are 
defined on large scales. Ultimately all flows appear to drain towards Shapley, 
although flows through the Perseus—Pisces filament take a circuitous route 
through the poorly studied Lepus region. 
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Extended Data Figure 2 | A representation of structure and flows due to 
mass within 6,000 km s™* (~80 Mpc). Surfaces of red and blue respectively 
represent outer contours of clusters and filaments as defined by the local 
eigenvalues of the velocity shear tensor determined from the Wiener filter 
analysis. Flow threads originating in our basin of attraction that terminate near 


the Norma cluster are in black, and adjacent flow threads that terminate at 
the relative attractor near the Perseus cluster are in red. The Arch and extended 
Antlia Wall structures bridge between the two attraction basins. The arch is part 
of a wall surrounding the Local Void. 
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Extended Data Figure 3 | One of three orthogonal views that illustrate 

the limits of the Laniakea supercluster. This SGX-SGY view at SGZ = 0 
extends the scene shown in Fig. 2 with the addition of dark blue flow lines away 
from the Laniakea local basin of attraction, and also includes a dark swath at 


SGY = 0 showing the region obscured by the plane of the Milky Way. As in 
Fig 2 in the main text, the orange contour encloses the inflowing streams, 
hence, defines the limits of the Laniakea supercluster containing the mass of 
10’” Suns and 100,000 large galaxies. 
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Extended Data Figure 4 | The second orthogonal view that illustrates the SGX = —4,750 km s`: the dark swath indicates the region obscured by the 
limits of the Laniakea supercluster. This is an SGY-SGZ slice at Milky Way. 
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Extended Data Figure 5 | The third orthogonal view that illustrates the SGY = +1,000kms_'. The Milky Way lies essentially in a plane parallel to the 
limits of the Laniakea supercluster. Here we show an SGX-SGZ slice at slice, at SGY = 0. 
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